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star's measurements of directed flow around midrapidity for tt*, K*, Kg, p and p in Au 
+ Au coUisions at ^snn = 200 GeV are presented. A negative v\ (y) slope is observed for most of 
produced particles (tt*, K^, Kg and p). In 5-30% central collisions a sizable difference is present 
between the v\{jj) slope of protons and antiprotons, with the former being consistent with zero 
within errors. The vi excitation function is presented. Comparisons to model calculations (RQMD, 
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UrQMD, AMPT, QGSM with parton recombination, and a hydrodynamics model with a tilted 
source) are made. For those models which have calculations of vi for both pions and protons, none 
of them can describe vi (y) for pions and protons simultaneously. The hydrodynamics model with 
a tilted source as currently implemented cannot explain the centrality dependence of the difference 
between the v\ (y) slopes of protons and antiprotons. 



PACS numbers: 25.75.Ld 

The BNL Relativistic Heavy Ion Collider (RHIC) 
was built to study a new form of matter known as 
the Quark Gluon Plasma (QGP) which existed 
in the universe shortly after the Big-Bang. At RHIC, 
two nuclei are collided at near light-speed, and the 
collision produces thousands of particles due to the 
significant energy deposited. The collective motion 
of the produced particles can be characterized Q by 
Fourier coefficients, 

Vn = {cosn{(j) - ip)) (1) 

where n denotes the harmonic, (/) and ip denote the 
azimuthal angle of an outgoing particle and reaction 
plane, respectively. The reaction plane is defined by 
the collision axis and the line connecting the cen- 
ters of two nuclei. Thus far, five of these coefficients 
have been measured and found to be non-zero at 
RHIC 0. They are directed flow wi, elliptic flow 
V2, triangular flow v^, the 4*'' order harmonic flow 
V4 and the 6*'' order harmonic flow wg. This paper 
will focus on the directed flow, the flrst Fourier co- 
eflieient. 

Directed flow describes the sideward motion of 
produced particles in ultra-relativistic nuclear col- 
lisions. It is believed to be generated during the 
nuclear passage time before the thermalization hap- 
pens, thus it carries early information from the col- 
lision 0-0] ■ The shape of directed flow at midra- 
pidity may be modified by the collective expansion 
and reveal a signature of a possible phase transition 
from normal nuclear matter to a QGP [sl-fioj. It 
is argued that directed flow, as an odd function of 
rapidity (y), may exhibit a small slope (flatness) at 
midrapidity due to a strong expansion of the fire- 
ball being tilted away from the collision axis. Such 
tilted expansion gives rise to anti-flow [1| or a 3'"'' 
flow component (not the third flow harmonic). 
The anti-flow (or the 3'"'' flow component) is perpen- 
dicular to the source surface, and is in the opposite 
direction to the bouncing-off motion of nucleons. If 
the tilted expansion is strong enough, it can even 
overcome the bouncing-off motion and results in a 
negative vi (y) slope at midrapidity, potentially pro- 
ducing a wiggle-like structure in vi{y). Note that 
although calculations Q for both anti-flow and 
3'"'' flow component are made for collisions at SPS 
energies where the first order phase transition to a 
QGP is believed to be the most relevant [l^, the di- 
rect cause of the negative slope is the strong, tilted 



expansion, which is also important at RHIC's top 
energies. Indeed hydrodynamic calculations for 
An -I- Au collisions at ^snn = 200 GeV with a tilted 
source as the initial condition can give a similar neg- 
ative vi{y) slope as that found in data. A wiggle 
structure is also seen in the Relativistic Quantum 
Molecular Dynamics (RQMD) model [l^l, and it is 
attributed to baryon stopping together with a pos- 
itive space-momentum correlation. In this picture, 
no phase transition is needed, and pions and nucle- 
ons flow in opposite directions. To distinguish be- 
tween baryon stopping and anti-flow, it is desirable 
to measure the vi [y) for identified particles and com- 
pare the sign of their slopes at midrapidity. In par- 
ticular, the observation of a centrality dependence of 
proton f 1 (y) may reveal the character of a possible 
first order phase transition [l^. It is expected that 
in very peripheral collisions the bouncing-off motion 
dominates over the entire rapidity range, and pro- 
tons at midrapidity flow in the same direction as 
spectators. In mid-central collisions, if there is a 
phase transition, the proton wi(y) slope at midra- 
pidity may change sign and become negative. Even- 
tually the slope diminishes in central collisions due 
to the symmetry of the collisions. 

At low energies, the E895 collaboration has shown 
that has a negative vi {y) slope around mi- 
dapidit y [l3| . while A and protons have positive 
slopes |14{ . This is explained by a repulsive kaon- 
nucleon potential and an attractive A-nucleon poten- 
tial. The NA49 collaboration [l^ has measured vi 
for pions and protons, and a negative vi(y) slope is 
observed by the standard event plane method. The 
three-particle correlation method t'i{3} p^ . which 
is believed to be less sensitive to non-flow effects, 
gives a negative slope too, but with a larger sta- 
tistical error. The non-flow effects are correlations 
among particles that are not related to the reaction 
plane, including the quantum Hanbury Brown- Twiss 
correlation , resonance decays [l^l and so on. At 
top RHIC energies, vi has been studied mostly for 
charged particles by both the STAR and the PHO- 
BOS collaborations [l9l - [2^ . It is found that vi in 
the forward region follows the limiting fragmenta- 
tion hypothesis [23| . and vi as a function of pseu- 
dorapidity (ry) depends only on the incident energy, 
but not on the size of the colliding system at a given 
centrality. Such system size independence of vi can 
be explained by the hydrodynamic calculation with 



a tilted initial condition The systematic study 



of vi for identified particles at RHIC did not begin 
until recently because it is more challenging for two 
reasons: 1) vi for some identified particles (for exam- 
ple, protons) is much smaller than that of all charged 
particles, 2) more statistics are needed to determine 
vi for identified particles other than pions. 

54 million events from Au + Au collisions at 
^snn = 200 GeV have been used in this study, all 
taken by a minimum-bias trigger with the STAR de- 
tector during RHIC's seventh run in year 2007. The 
main trigger detector used is the Vertex Position De- 
tector (VPD) [iJl. The centrality definition of an 
event was based on the number of charged tracks in 
the Time Projection Chamber (TPC) |25| with track 
quality cuts: I77I < 0.5, a Distance of Closest Ap- 
proach (DCA) to the vertex less than 3 cm, and 10 or 
more fit points. In the analysis, events are required 
to have the vertex z within 30 cm from the center of 
the TPC, and additional weight is assigned to each 
event in the analysis accounting for the non-uniform 
VPD trigger efficiency in the vertex z direction for 
different centrality classes. The event plane angle 
is determined from the sideward deflection of spec- 
tator neutrons measured by STAR's Shower Maxi- 
mum Detector inside the Zero Degree Calorimeters 
(ZDC-SMDs). Such sideward deflection of spectator 
neutrons is expected to happen in the reaction plane 
rather than participant plane, since the ZDC-SMDs 
are located close to beam rapidity. Being 6 units in 
r] away from midrapidity, ZDC-SMDs also allow a 
measurement of vi with minimal contribution from 
non-fiow correlations. The description of measuring 
vi using the ZDC-SMDs event plane can be found 
in [2l[. Particle Identification (PID) of charged par- 
ticles is achieved by measuring ionization energy loss 
(dE/dx) inside STAR's TPC, together with the mea- 
surement of the momentum (p) via TPC tracking. 
Track quality cuts are the same as used in [2^. In 
addition, the transverse momentum pT for protons 
is required to be larger than 400 MeV/c, and DCA 
is required to be less than 1 cm in order to avoid in- 
cluding background protons which are from knock- 
out/nuclear interactions of pions with inner detector 
material. The same cuts are applied to antiprotons 
as well to ensure a fair comparison with protons. 
The high-end of the pT cut is 1 GeV/c where protons 
and pions have the same energy loss in the TPC and 
thus become indistinguishable. For pions and kaons, 
PT range is 0.15 - 0.75 GeV/c and 0.2 - 0.6 GeV/c, 
respectively. Kg (—J- tt+tt") are topologically recon- 
structed by their charged daughter tracks inside the 
TPC [13. 

Results presented in the following figures contain 
only statistical errors. Results for pions, protons and 
antiprotons are not corrected for the feeddown from 
weak decay particles. The major systematic error in 
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FIG. 1: t;i for 7r±, K±, (left panel), p and p (right 
panel) as a function of rapidity for 10-70% Au + Au col- 
lisions at y^SNN = 200 GeV. The lines present the linear 
fit to the TT^, K^, Kg, p and p's vi{y) respectively. Data 
points around y = 0.29 are slightly shifted horizontally 
to avoid overlapping. 



determining the slope of vi (y) for identified particles 
is from the particle misidentification, which was eval- 
uated by varying the dE/dx cut. Another systematic 
error comes from the non-uniform pT acceptance, as 
vi{y) is obtained by integrating vi over the pr ac- 
ceptance which itself depends on the rapidity. This 
effect is non-negligible for protons and antiprotons 
at large rapidity. It is estimated by taking the dif- 
ference between slopes fitted with points integrated 
with PT acceptance at midrapidity and at large ra- 
pidity. In addition, some of the observed protons 
have originated from interactions between the pro- 
duced particles and the detector material, and such 
effect has also been taken into consideration. The 
total systematic uncertainty is obtained by adding 
uncertainties mentioned above in quadrature. There 
are also common systematic errors that should be 
applied to all particles: the uncertainty due to the 
first order event plane determination, which was es- 
timated to be ~ 10% (relative error) [2l|, and the 
uncertainty due to centrality selection, which was 
estimated to be ^ 4% (relative error) by comparing 
our charged ^1(77) slope to that from the RHIC run 
in 2004. Other systematic errors have been evalu- 
ated to be negligible. 

In Fig. [U vi{y) of tt^, K^, K^, p, and p are 
presented for centrality 10-70%. Following conven- 
tion, the sign of spectator vi in the forward re- 
gion is chosen to be positive, to which the mea- 
sured sign of vi for particles of interest is only rel- 
ative. Fitting with a linear function, the slopes are 
-0.15 ± 0.05(stat) ± 0.08(sys)(%) for the protons, 
-0.46 ± 0.06(stat) ± 0.04(sys)(%) for the antipro- 
tons, -0.27 ± O.Ol(stat) ± 0.01(sys)(%) for the pi- 
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FIG. 2: Model calculations of pion (left panel) and pro- 
ton (right panel) vi{y) for 10-70% Au -I- Au collisions 
at -y/sNN = 200 GeV. QGSM* model presents the ba- 
sic Quark- Gluon String model with parton recombina- 
tion [3y^. Hydro* model presents the hydrodynamic ex- 
pansion from a tilted source [lH . 



ons, -0.02±0.11(stat)±0.04(sys)(%) for the kaons 
and -0.17 ± 0.02(stat) ± 0.04(sys)(%) for the K%. 
The relative 10% common systematic error for all 
particles is not listed here. The vi{y) slope for the 
produced particle types {tt^, K^, and p) are 
mostly found to be negative at mid-rapidity, which 
is consistent with the anti-flow picture. In particu- 
lar, kaons are less sensitive to shadowing effects due 
to the small kaon-nucleon cross section, yet it shows 
a negative slope. This is again consistent with the 
anti-flow picture. Interestingly, vi (y) for protons ex- 
hibits a clearly flatter shape than that for antipro- 
tons. While mass may contribute to the difference 
in slope between pions and protons/antiprotons, it 
cannot explain the difference in slope observed for 
antiprotons and protons. Indeed, the observed vi for 
protons is a convolution of directed flow of produced 
protons with that of transported protons (from the 
original projectile and target nuclei), so the flatness 
of inclusive proton vi{y) around midrapidity could 
be explained by the negative flow of produced pro- 
tons being compensated by the positive flow of pro- 
tons transported from spectator rapidity, as a fea- 
ture expected in the anti-flow picture. 

In Fig. [21 pion and proton vi{y) are plot- 
ted together with five model calculations, namely, 
RQMD [13, UrQMD (pi, AMPT [H, QGSM with 
parton recombination [30| . and slopes from an ideal 
hydrodynamic calculation with a tilted source [Tl| . 
The model calculations are performed in the same pt 
acceptance and ccntrality as the data. The RQMD 
and AMPT model calculations predict the wrong 
sign and wrong magnitude of pion vi{y), respec- 
tively, while the RQMD and the UrQMD predict the 
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FIG. 3: Charged pions (solid stars), protons (solid cir- 
cles) and antiprotons (solid squares) ui (y) slope {dvi jdy) 
at midrapidity as a function of centrality for Au -I- Au 
collisions at a/snn = 200 GeV. 



wrong magnitude of proton vi{y). For models other 
than QGSM which has the calculation only for pi- 
ons, none of them can describe Vi (y) for pions and 
protons simultaneously. 

In Fig. |3l the slope of viijj) at midrapidity is pre- 
sented as a function of centrality for protons, an- 
tiprotons, and charged pions. In general, the magni- 
tude of the vi (y) slope converge to zero as expected 
for most central collisions. Proton and antiproton 
vi{y) slope are more or less consistent in 30-80% 
centrality range but, diverge in 5-30% centrality. In 
addition, two observations are noteworthy: i) the 
hydrodynamic model with tilted source (which is a 
characteristic of anti-flow) as currently implemented 
docs not predict the difference in vi (y) between par- 
ticle species [3l|. ii) If the difference between vi of 
protons and antiprotons is caused by anti-flow alone, 
then such difference is expected to be accompanied 
by strongly negative vi slopes. In data, the large 
difference between proton and antiproton vi slopes 
is seen in the 5-30% centrality range, while strongly 
negative vi slopes are found for protons, antipro- 
tons and charged pions in a different centrality range 
(30-80%). Both observations suggest that additional 
mechanisms than that assumed in [lll . [3ll | are needed 
to explain the centrality dependence of the difference 
between the vi{y) slopes of protons and antiprotons. 

The excitation function of proton vi{y') slope F 
{— dvi/dy' at midrapidity) is presented in Fig |4l 
Values for F are extracted via a polynomial fit of the 
form Fy' + Cy'^ , where y = y/ybeam for which spec- 
tators are normalized at ±1. The proton wi(y') slope 
decreases rapidly with increasing energy, reaching 
zero around ^snn = 9 GeV. Its sign changes to 
negative as shown by the data point at a/5 nn = 17 
GeV, measured by the NA49 experiment . A sim- 
ilar trend has been observed at low energies with a 
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FIG. 4: Proton vi{y ) slope {dvi/dy ) at midrapidity 
as a function of center of mass collision energy, where 

y = y/ybea,n- 



slightly different quantity d{px)/dy' [32|, [SJ]. The 
energy dependence of vi{y') slope for protons is 
driven by two factors, i) the increase in the number 
of produced protons over transported protons with 
increasing energy, and ii) the vi of both produced 
and transported protons at different energies. The 
negative vi{y') slope for protons around midrapid- 
ity at SPS energies cannot be ex plai ned by trans por t 
model calculations like UrQMD jUl and AMPT |29t. 
but is predicted by hydro calculations [1, Q. The 
present data indicate that the proton vi slope re- 
mains close to zero at -y/sNN = 200 GeV as observed 
at ^snn = 9 GcV and -^snn = 17 GeV heavy ion 
collisions. Our measurement offers a unique check 
of the validity of a tilted expansion at RHIC top 
energy. 

In summary, STAR's measurements of directed 
flow of pions, kaons, protons, and antiprotons for 
Au + Au collisions at y^snn = 200 GeV are pre- 
sented. In the range of 10-70% central collisions, 
vi{y) slopes of pions, kaons (K^), and antiprotons 



are found to be mostly negative at mid-rapidity. In 
5-30% central collisions a sizable difference is present 
between the vi{y) slope of protons and antiprotons, 
with the former being consistent with zero within 
errors. Comparison to models (RQMD, UrQMD, 
AMPT, QGSM with parton recombination, and hy- 
drodynamics with a tilted source) is made. Putting 
aside the QGSM model which has the calculation 
only for pions, none of the other models explored 
can describe vi{y) for pions and protons simultane- 
ously. Additional mechanisms than that assumed in 
the hydrodynamic model with a tilted source [HI ISH 
are needed to explain the centrality dependence of 
the difference between the vi{y) slopes of protons 
and antiprotons. Our measurement indicates that 
the proton's vi{y) slope remains close to zero for Au 
-I- Au collisions at -^/snn = 200 GcV. These new 
measurements on the particle species and centrality 
dependence of vi (y) provides a check for the validity 
of a tilted expansion at RHIC top energy. 
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